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Abstract 

We describe the analytic continuation of two-loop four-point functions with one off-shell external leg 
and internal massless propagators from the Euclidean region of space-like 1^3 decay to Minkowskian 
regions relevant to all 1 ^ 3 and 2 —> 2 reactions with one space-like or time-like off-shell external leg. 
Our results can be used to derive two-loop master integrals and unrenormalized matrix elements for 
hadronic vector-boson-plus-jet production and deep inelastic two-plus-one-jet production, from results 
previously obtained for three-jet production in electron-positron annihilation. 



1 Introduction 



In recent years, considerable progress has been made towards the extension of QCD calculations of jet 
observables towards the next-to-next-to-leading order (NNLO) in perturbation theory. One of the main 
ingredients in such calculations arc the two-loop virtual corrections to the multi leg matrix elements relevant 
to jet physics, which describe either 1 — > 3 decay or 2 ^ 2 scattering reactions: two-loop four-point functions 
with massless internal propagators and up to one off-shell external leg. 

Using dimensional regularization with d 7^ 4 dimensions as regulator for ultraviolet and infrared 
divergences, the large number of different integrals appearing in the two-loop Feynman amplitudes for 
2 — > 2 scattering or 1 ^ 3 decay processes can be reduced to a small number of master integrals. The 
techniques used in these reductions are integration- by-parts identities and Lorentz invariance 
computer algorithm for the automatic reduction of all two-loop four-point integrals was described in 

The use of these techniques allowed the calculation of two-loop QED and QCD corrections to many 
2^2 scattering processes with massless on-shell external particles [p[, which require master integrals 
corresponding to massless four-point functions with all legs on-shell ^,0. The results in B are given for 
all three physical Mandelstam channels, which are related by analytic continuation. In ||], the full set 
of two-loop four-point master integrals with one external leg off-shell was computed, for the kinematical 
situation of a 1 — > 3 decay, by solving the differential equations in external invariants Q fulfilled by these 
master integrals. These integrals were employed in the calculation of the two-loop QCD corrections to the 
e~^e~ 3 jets matrix element and to the corresponding helicity amplitudes in which can be expressed 
as a linear combination (with rational coefficients in the invariants and the space-time dimension d) of 
the corresponding master integrals. The 2^2 scattering processes related to e+e" — > 3 jets by analytic 
continuation and crossing are both of high phenomenological importance: hadronic vector-boson-plus-jet 
production and deep inelastic two-plus-one-jet production. 

The two-loop four-point functions with all legs on-shell can be expressed in terms of Nielsen's poly- 
logarithms [0, which have a well-defined analytic continuation The continuation of the master 
integrals in the on-shell case was discussed in detail in 0. In contrast, the closed analytic expressions 
for two-loop four-point functions with one leg off-shell contain two new classes of functions: harmonic 
polylogarithms (HPLs) and two-dimensional harmonic polylogarithms (2dIIPLs) Accurate nu- 
merical implementations for HPLs and 2dHPLs |14| are available. The implementations apply to HPLs 



of arbitrary real arguments 13 1, but only for a limited range of arguments of the 2dHPLs. This range 



of arguments corresponds precisely to the Euclidean region (space-like 1 — > 3 decay) for the corresponding 
ofF-sliell master integrals ^ . Continuation from the Euclidean region to any physical (Minkowskian) region 
requires in general the analytic continuation of the 2dHPLs outside the range of arguments considered in 
[0. Only in the special case of a time-like 1^3 decay (relevant to e+e~ 3 jets), which corresponds 
to the simultaneous continuation of all three external invariants from Euclidean to Minkowskian values, 
this continuation can be carried out by simply replacing an overall scaling factor, while preserving all 2dH- 
PLs (which depend only on dimensionless ratios of the invariants). These results were used in |^] for the 
calculation of the two-loop matrix elements for e^e^ ~> 3 jets. 

It is the aim of the present paper to derive the relations for HPLs and 2dHPLs needed for the analytic 
continuation of the two-loop four-point master integrals of Q to the kinematics of all the 2^2 scattering 
reactions with one off-shell external leg, working out real and imaginary parts explicitly. We also provide 
the algorithms for expressing them in terms of HPLs and 2dHPLs whose arguments lie within the range 
covered by the numerical routine of 1 14 1 , which allows them to be evaluated numerically. Given that the 
matrix elements are linear combinations of the master integrals, this will in turn allow us to determine the 
matrix element for all 2 ^ 2 reactions related to e+e~ 3 jets by crossing. 

In the context of the next-to-leading order (NLO) corrections to jet observables, the first calculation 
was also for the kinematics of the 1 — > 3 decay, relevant to e^e~ — > 3 jets. The one-loop matrix element 
obtained in this calculation contained only logarithms and dilogarithms, which have a known analytic 
continuation. The results of |15| were continued to the kinematic situation relevant to hadronic vector- 
boson-plus-jet production in |16|, and to deep inelastic two-plus-one-jet production in jl^,|l^. The analytic 
continuation procedure used in these calculations is documented in detail in . In particular it is already 
observed at the one-loop level that the kinematic region relevant to the hadronic vector-boson-plus-jet 



1 



production is free of kinematic cuts, while two kinematic cuts are present inside the region of deep inelastic 
two-plus-one-jet production. We shall see below that the same pattern is preserved at the two-loop level. 

The paper is structured as follows. In Section we define the kinematical variables used to describe 
two-loop four-point functions with one off-shell leg and describe the ranges of variables relevant to each 
different process by making a detailed decomposition of the kinematic plane. The basic features of analytic 
continuation are discussed in Section ^ on the example of the continuation of the 1-^3 decay from the 
Euclidean to the Minkowskian region. In Sections ^ and |[ we derive the algorithms for the analytic 
continuation in one invariant for a time-like and a space-like off-shell leg. Section ^ contains conclusions 
and an outlook. Finally, the Appendix recalls definitions and main properties of the HPLs and 2dHPLs. 

2 Kinematic regions and notation 

We label the external momenta on the four-point functions with pi, p2, Ps and p4, and write momentum 
conservation a,s pi + p2 + Ps + P4 — , with on-shell conditions pf = for i = 1,2,3 while p4 is off-shell 
Pi = (^Pi ^ P2- Ps)^ = 9^ 7^ 0. We further define 

Sl2 = {Pl +P2)^ , Sl3 = {Pl +P3)^ , S23 = {P2 + Psf , (2.1) 

SO that the Mandelstam relation reads 

si2 + si3 + S23 = • (2.2) 

We will use the metric in which time-like invariants are positive. The kinematic plane defined by S12, S13 
and S23 is shown in Fig. ^ where equilateral (non-Cartesian) coordinates were used to display the symmetry 
in the three invariants. The lines indicate the locations of potential cuts in the four-point functions. For 
later use, all regions defined by these cuts are labelled as (la), (lb), . . ., (4d). 

The regions physically relevant to different processes are displayed in Fig. ||. In e+e" 3 jet (3j) 
production, is time-like (hence positive) and all Sij are positive as well. The relevant region is the inner 
triangle of the kinematic plane, as shown in Fig. |^(a). This inner triangle corresponds to region (la) in 
Fig. |l|. To indicate space-like {q^ < 0) and time-like {q^ > 0) kinematics for a region under discussion, 
we will use in the following the subscripts "— " (space-like) and (time-like). The region relevant to 
3j-production is thus denoted by (la)+. 




s,, = 



Figure 1: Kinematic plane in terms of Lorentz invariants = (pi+Pj)'^ displayed in equilateral coordinates. 
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(a) s„ = 





Figure 2: Regions of the kinematic plane relevant to (a) e'^e 3 jets {q^ > 0), (b) {V + 1) jet production 
{q^ > 0) and (c) deep inelastic (2 + 1) jet production {q^ < 0). 

Vector-boson-plus-jet {V + Ij) production at hadron colliders and deep-inelastic two-plus-one-jet (DIS- 
(2 -|- production are described by three subprocesses each (corresponding to the si2, Si3, S23 channels, 
which all contribute to these final states). 



For V + Ij production, q is time-like, and for pi + pj Pk + Pi the invariants fulfil 

q^ > , Sij > q^ > Q , s-jk < , Sjfc < , 



(2.3) 



where {i,j,k) stand for the three non-ordered permutations of (1,2,3). The relevant regions are shown in 
Fig. ||(b) and correspond to the regions (2a,3a,4a)+ of the kinematic plane. Fig. |^. 

Finally, for DIS-(2 + l)j production, q^ is space-like (hence negative) and for P4 + Pk Pi + Pj the 
invariants fulfil 

q^ <0 , Sy > > , sjk < 0, s,k < , (2.4) 

where {i,j,k) stand again for the three non-ordered permutations of (1,2,3). We display the relevant 
kinematic regions in Fig. ^(c). Each of those regions cannot be identified with a single region in the 
kinematic plane. Fig. |l], but is instead patched together from four regions, (Id, 2c, 3b, 4d), (lb, 2d, 3c, 4c) and 
(lc,2b,3d,4b). 

It is customary to introduce the dimensionless variables 



S12 



Sl3 



£23 
«2 



(2.5) 
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Figure 3: Kinematic plane in Cartesian coordinates. 



for which (2.2) reads 

X + y + z = l. (2.6) 

As independent variables, we will use mainly q^, y and z, with x therefore given as x = 1 — y — z. We 
will further represent the various kinematical configurations, at given q^, in the Cartesian y, z plane. For 
convenience of later use, we represent in Fig. ^ the whole y, z plane properly partitioned in all the regions 
which will be of interest later. The labelling of the regions is as in Fig. |^. 

When the three regions (3j, V + Ij, DIS-(2 + l)j), are superimposed, regardless of g^, they cover the 
entire (y, z) plane. 

Despite using mainly g^, y and z as independent variables, all analytic continuations are to be carried 
out in the original Mandelstam variables Si2, Si3 and S23, adding to them, when time- like, the usual 
infinitesimal imaginary part with positive sign, Sij + ie. Indeed, no definite sign can be attributed a 'priori 
to the imaginary parts of the dimensionless invariants a;, y and z, given the constraint x + y + z = 1. In 
practice, this means that any function of y and z that is to be continued analytically, has to be expressed 
first in terms of S12, S13 and S23, then continued taking correct account of the well defined imaginary 
parts of the s,j , then finally re-expressed in and y, z or other dimensionless variables appropriate to the 
considered region. 



3 Basics of analytic continuation 

The problem of analytic continuation is best discussed starting from the unphysical Euclidean case in which 
<Z^, S12, S13 and S23 are all space-like (hence negative). Let us give them the values 

with 

Q2 > , cr„- > , (T12 -I- ai3 + 0-23 = ■ (3.2) 
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All master integrals are indeed real in this configuration. According to Eq. (2.5), we introduce 



_ £12 _ CT12 _ fl3 _ 0^3 _ £23 _ O23 

Q-^ q-^ (4^ q^ 



and observe that, because of Eq. (3_^), we have (0 < y < 1, < z < 1 — y) or equivalently (0 < z < 1, < 
y < 1 — z), i.e. (y, z) are in the region (la) of Fig. ^. Any of the master integrals evaluated in say 
^(si2, si3, S23), can be written as 



*(S12,S13,S23) = ^QT^\^^-^ 
2\Oi^ ( ^13 ^23 



(-g^)"$ 

(-g2)"cO(y,z) , (3.4) 



where use is made of Eq. ( |3.3D and a is an exponent accounting for the mass dimension of the integral: 
a — d — r + s, where d is the continuous space-time dimension used in dimensional regularization, while 
the integers r and s denote the number of propagators and scalar products in the integral. For a two-loop 
seven-propagator master integral (e.g. a scalar double box or crossed box integral), with all propagators 
raised to unit power and no scalar products, one thus has a — d—1 — —3 4" (d — 4). The function $(y, z) is 
in general a Laurent series in (d — 4), whose coefficients are combinations of simple algebraic factors in y, z 
times HPLs of various weight and argument z and 2dHPLs of various weight and argument y, with indices 
depending on z. 

Note that the dimensionless ratios y, z are both real and positive, and lie within the boundaries < 
z < 1, < y < 1 — z. In that range of values of the arguments, all the HPLs and 2dHPLs are analytic and 
real - as expected, of course, as all the kinematical variables are Euclidean. This region coincides with the 
triangle (la)_ of Fig. |. 

As a consequence of the symmetry and analyticity properties of the HPLs and 2dHPLs, within the 
region (la)_ we can exchange the roles of y and z and re-express the master integrals as 

*(si2,si3,S23) = (-<z')"$'(z,y) , (3.5) 

where $'(z,y) consists of HPLs of argument y and 2dHPLs of argument z with indices depending on y, 
all within their analyticity region. The transformation from $(y, z) to $'(z, y) can be implemented, within 
the triangle (la) of Fig. ||, by re-expressing the HPLs H(m; z) and 2dHPLs G(m(z);y) as a combination 
of HPLs H(m; y) and 2dHPLs G(m(y); z), by means of the 'interchange-of-arguments' procedure described 
in Appendix A. 2. 2 of [||. All algorithms used for this procedure (and for all subsequent transformations 
derived in this paper) were coded in FORM |l^ . 

Similarly, within the region (la) we can replace z = l — x — ybyx, obtaining another representation of 
the master integrals 

*(si2, si3, S23) = (-g2)"$"(y, x) , (3.6) 

where $"(y,a;) consists of HPLs of argument x and 2dHPLs of argument y, with indices depending on x, 
all within their analyticity region. The above transformation, which amounts to re-expressing the HPLs 
H(to; z) and 2dHPLs G(TO(z);y) appearing in the original <i>(y, z) in terms of HPLs H(n;x) and 2dHPLs 
Q{n(x)\ y), with cc = 1 — y — z, can be implemented as a result of the combination of the 'interchange-of- 
arguments'-procedure described in Appendix A. 2. 2 of |^ with the reflection algorithm derived in Section 5 
of 0. 

The notation for the 2dHPLs used in the present work is shortly recalled in the Appendix. It is the 
same as was introduced in [Q , which is however different from the notation originally proposed and used 
in j|] to present the two-loop four-point master integrals. The notation of was already employed in [|j 
to represent the result for the two-loop QCD corrections to the 3j matrix element. Detailed transformation 
rules between the two notations can be found in . A summary of definitions and properties of HPLs 
and 2dHPLs is provided in the appendix of [b|. 
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In the 3j case, all the kinematical variables are time- like, hence positive; the proper analytic continuation 
is obtained by starting from the Euclidean case and giving to each positive variable (he nce t o all the variables 
in this case) a small positive imaginary part (the usual +ie prescription). Equation (3.4) then becomes 



*(si2 + ie, si3 + if-, S23 + if) = (-9^ - «e)"$ 



The real parts of (si3 + ie)/{q^ + if) and (523 + if)l{q^ + if)-, 



si3 + »e 523 + 



Sl3 



S23 



(3.7) 



(3.8) 



lie in the region 0<y<l,0<z<l — y, the triangle (la) of Fig. |^, where all the HPLs an d 2d HPLs are 
analytic and real. The e ^ limit is then trivial (i.e. the it can be simply ignored) and Eq. (3.7) becomes 



*(si2 + it, Si3 + it, S23 + if) = {-q^ - ifT^iv-, z) , 



(3.9) 



where $(?/, z) is exactly the same as in the fully Euclidean case. The continuation from (la)_ to (la)+ 
performed here leaves $(2/, z) unchanged. The imaginary parts of the 3j master integrals (which are of 
course complex) are entirely due to the factor {—q^ — it)". The integer part of a does not matter, while 
the expansion in {d — 4) gives 



i-q'-^f) 



d-A 



= 1 + (d - 4)(ln(g2) - + i(d _ 4)2(ln(g2) - i^f + .. 



(3.10) 



It is worth recalling that the master integrals, as well as the physical matrix elements, develop polar 
singularities around d = 4 in dimensional regularization, so that these higher-order terms in the {d — 4)- 
expansion and their imaginary parts become of actual importance. 

It is to be emphasised here that in H], strictly speaking, the master integrals ^'(si2, S13, S23) were never 
directly evaluated in the fully Euclidean region; but as_$(y, z) is the same functi on, and for the same range 
of arguments, in both the fully Euclidean case, Eq. ( |3.4| ), and the ij case, Eq. (3^), we can define in that 



way (i.e. by just replacing in the overall scale factor {—q —it)"" by (Q )") the Euclidean master integrals in 
terms of the master integrals given in ||^ for the 3j case. From now on, we will therefore take the Euclidean 
master integrals as known, and will show how to get by analytic continuation the master integrals in the 
various kinematical regions of physical interest. 



4 Analytic continuation in one invariant for time-like 

To continue from region (la)_ to regions (2a,3a,4a)+, which are relevant to 4- Ij production at hadron 
colliders (with a time-like momentum of the vector boson V) , it is necessary to continue simultaneously in 
one of three Lorentz invariants s,;j and in the vector boson virtuality q^ . 
We start by discussing the continuation from (la)_ to (4a) + Q 

In region (4a)+ of Fig. |^, relevant for V +lj production in the momentum arrangement p2 +P3 ~* Pi +P4, 
(1 < z < 00, 1 — z < y < 0) or (—00 <y<0, 1 — y<z< 00) and 

q^ > 0, S12 < 0, si3 < 0, S23 = q^ - S12 - si3 > > , (4.1) 

so that for the analytic continuation q^ and S23 must be given infinitesimal imaginary parts +it. From 
Eq. (3^) one gets in this case 

«'(si2, Si3, S23 + ^f) = {-q^ - *e)"<f f ^f^) = - + ^ - (4-2) 

\q'^ + It q^ + It J 



^Notice that the analytic continuation from (la)_ to (3a)+ as outUnod in the appendix of [Q| is unnecessarily complicated. 
Moreover, Eq. (A. 30) in [M, which is relevant in this context, contains a misprint: all i's should read —i. 
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where, because of Eqs. (11) and of the very definition of y, z, Eqs. (2^), we have used 



o , . = ^ + »e = y + *e , n , . = ^ - «e = z - . (4.3) 

We introduce new dimensionless variables u, v through the relations 

y q2 1 

M4a = W = = , Via=V = = - , (4.4) 

8-23 Z S23 Z 

so that 

y = , z = - . 4.5 

V V 

As y and z span the region (4a) of Fig. |3[ we find that u, v vary in the ranges < f < 1 and < u < I — v, 
i.e. the above parametrisation maps region (4a) into region (la). 
Equation (4.2) then reads 

*(si2, si3, S23 + ie) = {-q^ - ie)"$ ( -- +ie,- - ie] . (4.6) 

\ V V J 

Given the above prescription, it is relatively straightforward to determine the proper analytic continu- 
ation of the various HPLs and 2dHPLs appearing in the analytic expression $(y, z) of any master integral 
and then to express them as HPLs of argument v and 2dHPLs of argument u and indices depending on 
V, with {u,v), as already observed, within the analyticity triangle of the functions. Let us start from the 
HPLs and 2dHPLs of weight w equal to 1, which are just logarithms. According to the definitions of the 
appendix, one has 

H(0;z — ie) = log ^ ie^ = —\ogv 

= -H(0;«), 

H(l; z — ie) = — log ^1 — - + ie^ = — log(l — v) + \ogv — in 

= H(l;i;) +H(0;i)) -iTT . (4.7) 

Notice that for H(0; z), the ze-term can safely be ignored. 

The 2dHPLs at w — 1, again according to the definitions of the appendix, are continued as 

/ u \ / u 

G{0;y + te) = log + ie = log - 

\ V / \v 

= G(0;u) -H(0;u) +i7r 

(11 
1 + -- 
V 

= G{-v;u) 

/l + --ie-i , , 

G{1 -z + ie;y + ie) = log ^ ^ = log 

V 1 - + «e 



= 0(1 -v;u) 

Gi-z + ie;y + ie) = log (^^^^^^^^^ j = log(l - u) 

= G{l;u). (4.8) 

It should be noted that no definite imaginary parts can be assigned to the arguments of the logarithms in 
G(l — z + ie;y + ie) and G(— z + ie;y + ie); but as the arguments in both cases remain positive and within 
the analyticity region of the functions, the result is anyhow well defined. 

Using the above formulae, all imaginary parts of the higher weight HPLs and 2dHPLs are fixed, since 
these functions can be derived iteratively, starting from the w = 1 functions, as will now be shown. 
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If, in an HPL of weight w higher than 1, aU the indices are equal to 1, that HPL is just equal to 1/wl 
times the w-th power of H(l; z — ie), already seen. When the indices a re no t all equal to 1, one can first 
separate all leftmost (l)'s in the index vector, using the product algebra ( A. 9 ), such that the leftmost index 
is always a (0). Any HPL can then be written as [|3| 

H(0,M-ie) = H(0,&;l-ie)+ / dz' ^ll{0,b; z') 

= H(0,fe;l)+ / dz'— H(0,6;z'-ie) 
Ji dz' 

= H(0,&;l)-^"^H(^6;i-ze^ , (4.9) 

where in the last step the integration variable v' — 1/z' was introduced. The expression for H(fe; 1/v' — ie) 
in terms of H(c; v') and its proper imaginary part is of lower weight than H(0, b; z — ie) and thus already 
known in an iterative bottom up approach in the weight w. As an example of this transformation, one finds 

H(0, 1, 1; z - le) = H(0, 0, 1; v) - H(0, 1, 1; v) + H(0; v) - H(0, 1; v) ~ ^RiO; «)H(0; v)^ + (3 

-H(0, 1; v) ~ iH(0; t;)H(0; + y ) • (4.10) 

To perform the analytic continuation of the 2dHPL in y, one first separates off all rightmost (0)'s in the 
index vector by applying the product algebra. The remaining 2dHPL (the dependence of the indices a, h 
OM z — Xjv ~ ie is understood for short) can then be written as 

G(a,6; — + = / dy'— G(a, 6; y') 



V ^ Jo dy 



^ " dy'^G{a,b;y' + ie) 

dy'g{a,y' + ie)G{b-y' + ie) 

' du' -g (a, -- + ie) G (b- -- + ze ) , (4.11) 



V \ V 

where the expression for G(6; —u'/v-\-ie) in terms of G(c; u') and its proper imaginary part is again of lower 
weight and thus known in an iterative bottom up approach in w. In the rational fractions the ie does not 
matter and their expressions are 

-g(o;--) = -g{0;u'), 

V \ V J 

-g(i;-— ) = -gi-v;u') , 



-g l-z; = -g(l-i;;-u), 



V \ V 

1 / u' 

-g -z; 

V \ V 



l{l;u'), (4.12) 



such that the w'-integral in ( 4.11 ) yields a 2dHPL of argument u. 

From the above, it becomes immediately clear that only a 2dHPL with trailing (0)'s in the index vector 
acquire an imaginary part when continued from (la)_ to (4a)4_. Moreover, any 2dH PL from (la)_ without 
trailing (0)'s is identified with a single 2dHPL in (4a) + (since the w — 1 functions ( ^^ ) are identified on a 
one-to-one basis). 
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An example of the continuation of a 2dHPL is 



G(0, 1,1-0,0;?/ + ie) = -H(0, v) G(0, -v, l-v,u) + G(0, -v, l-v,0,u) + iir G(0, -v, l-v;u) . (4.13) 



To summarize, master integrals and matrix elements, which ar e gi ven in (la)_ in terms of HPLs H(m; z) 
and 2dHPLs G(m(z); y) can be continued to (4a)+, by using Eq. ( |4.5| ), where they are expressed in terms of 
HPLs H(m;ti) and 2dHPLs G(to(i;);m). Given the definitions u — ~y/z and v — 1/z, one finds < v < 1, 
0<u<l — win (4a)^uch that the above HPLs and 2dHPLs are real, and can be numerically evaluated 
with the routines of |l^,|l4|. Imaginary parts were made explicit in the analytic continuation. 

The other two momentum arrangements relevant to vector boson production at hadron colliders are 
Pi + P2 — > P3 +P4, corresponding to region (2a)-|_ and pi + ps P2 +P4, corresponding to region (3a)+. 

In the region (2a)+, we have {—oo < z < 1, — oo < y < 1) and 



> 0, Si2 = g - Sl3 - S23 > r > 0, Si3 < 0, S23 < , 

so that and Si2 must, for the analytic continuation, be given infinitesimal imaginary parts 



(4.14) 



In close analogy with the discussion for the region (4a) +, we start from Eq. (3.4), which now becomes 

Sl3 S23 



^-(512 + ie, Si3, S23) = i-q - ier<i> 



q"^ + ie' q"^ + ie 



= {-q - ie)°'^{y + ie,z + ie) 



with 



Sl3 



S13 

= — +ie^y + ie , 



S23 



q^ + le q- q^ + le 

we then introduce new dimensionless variables U2a7W2a as 



W2a 



so that 



£13 
S12 

y 



y 



i-y 

W2a 
W2a 



V2a. 



1 - M2a 



£23 
•2 



q ^ 

S12 

- V2a 



1-y- Z 



V2a. 



(4.15) 

(4.16) 

(4.17) 
(4.18) 



and M2a,W2a vary in the ranges < z)2a < 1 and < U2a < 1 — W2a, i-e. the above parametrization maps 
region (2a) into region (la), and Eq. (4.15) becomes 



*(si2, Si3, S23 + ie) = {-q^ 



)"$ 



U2a 
V2a 



1 - W2a 



W2a 



V2a. 



(4.19) 



The separation of the real and imaginary parts of the HPLs and 2dHPLs of the above arguments — M2a/w2a + 
ie, (1 — U2a — ti2a)/w2a + «£, and their expression in terms of HPLs and 2dHPLs of ar gume nts U2a, W2a and 
(i7r)'s can then be carried out by a suitable extension of the derivation of Eqs. (4.7)~(4.13). 

Alternatively, we can use the representation ( [3.(^) of the master integrals, i.e. we can first transform 
the expression for the Euclidean master integral ( |3.4| ) by rewriting $(y, z) in terms of a; = 1 — y — z, so 
obtaining the function ^"{y,x) defined as 



$"(y,x) ^ <i>{y,l~x-y) 



For the continuation to (2a)+, Eq. (3.6) reads 

Sl3 



*(si2 + ie, si3, S23) = ("9 - ie)"* 



S12 



with 



Sl3 



S13 

-^+ie = y + ie , 
q^ 



q^ + ie' q^ + ie 
S12 + 



{-q^ - ie)"$"(j/ + ie,x- ie). 



S12 



The dimensionless variables U2a7 W2a of Eq. (4.17) can also be written as 

^2 



U2a 



£13 
S12 



y_ 

X 



W2a 



S12 



1 

X 



(4.20) 

(4.21) 
(4.22) 

(4.23) 
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so that 



and Eq. ( 4.21 ) becomes 



U2a 
V2a. 



W2a 



*(si2 + le, si3, S23) = (-r - *e)"*" i€ 



U2a 



1 



(4.24) 



(4.25) 



The same analytic continuation formulae as applied above for the continuation of Eq. ( |4.6D to (4a) + can 
therefore be used in this case as well, if allowance is made for the formal replacement z x and v — s- V2a- 
In the region (3a)_|_ , corresponding to pi + P3 ^ P2 + Pi, (—00 < z < 0, I — z < y < 00) and 



> 0, S12 < , si3 = - S12 - S23 > 9^ > , S23 < , 
for the analytic continuation and S13 must be given imaginary parts +ie, 

*(si2, S13 + ie, S23) = {~q^ - ze)"* f ^f:^, = i-q' - i^THy -ie,z + ze) 

\ q^ + le q^ + le J 

We introduce new dimensionless variables U3a,t'3a as 



W3a 



£23 

Sl3 



si3 y 



y 



(4.26) 
(4.27) 

(4.28) 
(4.29) 



so that 

1 ^ _ _M3a 

l'3a ' V3a 

i*3a, i^3a Vary in the ranges < < 1 and < < 1 — W3a7 i-e. the above parametrization maps region 
(3a) into region (la), and the proper analytic continuation is given by 



*(si2, S13 + ie, S23) - (-9' - *£)"$ ie,- — 



1 



U3a 



W3a 



(4.30) 



In close analogy with the (2a) + case, the separation of the real and imaginary parts of the HPLs 
and 2dIIPLs of arguments l/w3a — le, — M3a/w3a + *e and their expression in terms of HPLs and 2dHPLs of 
arguments u^a., v^g. and (z7r)'s can then be carried out by a suitable extension of the derivation of Eqs. ( [4.7| )- 
( 4.13 ) previously established for the region (4a) +. 

Alternatively, we can start from the expression Eq. (3_^) for the Euclidean master integral obtained by 
crossing the arguments: 



«'(si2,si3,S23) = (-g^)"$(y,z) 
= (-g2)"$'(z,y) 



(4.31) 



The expression (4.31) then reads 

^"(312,313 + ie, S23) 



{-q^ - ie)°$' 
{~q^ -ie)"$'(z + ie,y-ie) 



523 gi3 + »e 



U3a 



= i-q' - ie^^' { - —+ie, ie 



V3a 



(4.32) 



with U3a,V3a givcn by Eq. ( 4.28 ). To this expression, we can again apply the same analytic continuation 
formulae as above, when allowance is made of the formal replacement y ^ z and u M2a, v V2a- 
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5 Analytic continuation in one invariant for space-like 



To continue from region (la)_ to regions (lbcd,2bcd,3bcd,4bcd)_ , which are relevant to deep inelastic two- 
plus-one-jet production, it is necessary to continue one of the three Lorentz invariants Sij to the time-like 
region, while not altering the negative sign of . 

In total, there are twelve regions in the kinematic plane that are relevant to DIS (2 -I- l)j-production. 
It turns out that the analytic continuation to all these regions can be obtained by deriving continuation 
formulae to four regions (which we take to be (ld)_, (4d)_, (4b) _ and (3c)-), while the continuation to the 
remaining eight regions is then obtained using crossings of arguments, as described in the previous section. 
We establish the continuation formulae for all these cases in this section. 



5.1 Continuation from (la) to (Id) and to (lb) , (Ic) 

In region (ld)_, which contributes to DIS in the momentum arrangement Pi+pi P2 + Ps (but does not 
cover the full phase space available for this reaction), we have 0<?;<1, Q > z > — y, i.e. — 1 < z < 
0, —z<y<l, and 



q-" < 0, si2 < 0, 



q < Si3 < 0, -Si3 > S23 > 



(5.1) 



such that only S23 needs to be assigned an infinitesimal imaginary part +ie. Owing to (5.1 

Sl3 



one has 



S12 + S13 + S23 + 



y + ie , 



Equation (3^) therefore reads 
^(si2,si3,S23 + 



(-g2)"$ 



S23 + »e 
S12 + .si3 + S23 + 



Sl3 



S23 + ?e 



(5.2) 



S12 + si3 + S23 + «e S12 + si3 + S23 + ie 



= i-q^r^{y + ie,z-ie) 



(5.3) 



Note that the sign of the imaginary part associated with y, as will be clear from the following, actually 
plays no role in the assignment of imaginary parts to the 2dIIPLs, since y remains in the range < y < 
1 < 1 — z in (ld)_, which is free from cuts. 

In region (ld)_, we introduce the new variables 



''Id 



S12 + S23 



£23 
«2 



(5.4) 



so that 



Sl3 , S23 

y = = 1 ^ : Z = 

q Sl2 + Sl3 + S23 



(5.5) 



in the region (ld)_, (r, s) fulfil 0<»'<l,0<s<l — r, orO<s<l,0<?'<l~s, thus mapping (Id)- 
onto (la)_. 



The generic expression for a master integral, Eq. (3.4), then reads 

*(si2, si3, S23 + ie) = (-q2)"$ {y + ie, z - ie) = {-q^)"<S> {1 - r + ie, -s - ie) . 
In terms of these variables, the HPLs of ui = 1 are continued as 

11(0; z — ie) = log (— s — ie) — log (s) — iir 

= H(0;s)-i7r 

11(1; z-«e) = - log (1 + s + ie) 

= -H(-l;s), 



(5.6) 



(5.7) 
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and the 2dHPLs at = 1 as 



G{0;y + ie) = log(l-r + ie) 

= G(l;r) 

G(l;y + ie) = log (r - ie) 

= G(0;r) 

G(l-z + ie;y + ie) = log (^^-^^^^^-t^^ = log (^^^^ + log(s) - log(l + s) 

= G(-s;r) +H(0;s) -H(-l;s) 

/ \ — T ~\~ Z€ \ ( 1 — T 

G(-z + ie;j/ + ie) = log 1 — — = log 1 

V s + le / \ s 

= G(l-s;r)-H(l;s)-H(0;.s) + i7r, (5. 



where we have used 



1 — r + ie 1— r 
: — = fe 



as < s < 1 — r in the region (ld)_. Using the above formulae, all imaginary parts of the higher weight HPL 
and 2dHPL are fixed, since these functions can be derived iteratively by integrating the w = \ functions. 

The HPLs with w > \ are obtained by first separating all rightmost (0)'s in the index vector using 
the product algebra [ p^ . The continuation to z — ie = — s — (flipping the sign of the argument) is then 
carried out (when the rightmost index a\ is different from 0) as described in ]l3t : 

H(a2,ai;2-ie) = (-l)"^+"^H(-a2, -ai; s) , 
H(a3,a2,ai;^-ie) - (-l)'^^+'^^+'^^'H(-a3, -aa, -ai; s) , 
H(a4,a3,a2,ai;^- ie) = (-l)'^^+''^+'^^'+'^''H(-a4, -aa, -03, -ai; s) . (5.9) 

Notice that the it in the argument is relevant only to the HPLs with rightmost (0)'s, which acquire an 
imaginary part in this transformation. 

To perform the continuation of the 2dHPLs with weight > 1, one can proceed by induction on 
the weight w, starting from the w — \ formulae (5.8). At w > 1 one first separates off all leftmost (1) 
components of the index vector by using the product algebra of the 2dHPLs, so that one has to consider 
only vector indices of the form (a(z),6(z)), where a(z) stands for one of the values (0, 1 — z, while 
h{z) can contain the index (1) as well. By using the very definition of the 2dHPLs, one can write in full 
generality 

Q(a{z),h{z)-y) = C dy'g{a{z);y')G{b{z);y') 



= G(a(z), b{z)- 1) + dy'g{a{z)-y')G(b{z)-y') 

= G{a{z),b{z);l)~ f ^ dr'g{a{z);l ~ r')G(b{z);l - r') , (5.10) 
"'0 

where the new integration variable r' = 1 — y' has been introduced. Let us recall that G(a(z),&(z); 1) is 
finite as the leftmost index a{z) is different from 1. 

The above formula is well suited for the continuation to (z = — s — ie, y = I — r + ie): 

G(a(— s — ie), b{—s — ie); I — r + ie) — 

= G{a{-s -ie),b{~s -ie);l) - dr'g{a{-s - ie);l ~ r')G{b{-s - ie);l - r') 



= G(a(-s - ie),b{~s - ie); I) - I dr'g{a{~s - ie): 1 - r' + ie)G(6(-s - ie): 1 - r' + ie) . (5.11) 
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The values at y — 1 can be evaluated by expressing G(a(z),6(z); 1), for < z < 1, in terms of HPLs of 
argu ment z, and then continuing the resulting expression to z ^ ie — ~s — ie as already discussed above, 
Eqs. (|5.9|) . For evaluating the integral, note that the integration variable runs in the region < r' < r < 
1 — s, so that the expressions of the rational fractions g(a(— s — ie); 1 — r' + ie) are 

g{0;l -r' + ie) = -g(l;r'), 

g(l;l-r' + ?:e) = -g(0;r'), 

g{l + s + ie;l-r' + ie) = -g(-s;r'), 

g{s + ie;l-r' + ie) = -g(l-s;r'), 

where the ie can be dropped as they are irrelevant in the considered region of r' . Finally, the last term 



appearing in ( 5.11 ), G(6(— s — ie): 1 — r' + ie) is of weight w — 1, so tha t its expression in terms of 2dHPLs 



depending on (s, r') is also known on that region, and the r'-intcgral in ( 5.11 ) can be immediately evaluated 
in terms of 2dHPLs of argument r. 

An example for the continuation of a 2dHPL is 

G(-z + ie,l^ z + ie;y + ie) = G(l - s, -s; r) + [H(0; s) - H(-l; s)] G(l -s;r) + H(-l, 0; s) 

+H(-1, 1; s) + H(l, -1; s) - H(l, 0; s) - ^ - H(-l; s) , (5.12) 



where use has been made of 



^2 



G(-z + ie,l-z + ie: 1) = H(-l, 0; s) + H(-l, 1; .s) + H(l, -1; s) - H(l, 0; s) iir H(-l; s) . (5.13) 

6 

This example also illustrates the major new feature encountered in the continuation from (la)_ to (ld)_ 
(or to any of the regions labelled by (b,c,d) in Fig. ||): the cut structure of the boundaries of the (b,c,d)-type 
regions does not reproduce the cut structure of the boundaries of the (a)-type regions. At each corner of 
any (a)-type region, three cuts intersect, while all (b,c,d)-type regions have one corner with only two cuts 
intersecting, plus other corners (if any) with three cuts intersecting. In the case of the region (ld)_, it is 
the lower right corner, which touches only two cuts {y — \ and y = — z), while the two upper corners touch 
three cuts each (z = 0, y = and y = —z for the upper left corner and z = 0, y = 1 and y = 1 — z for the 
upper right corner respectively). Moreover, the (a)-type regions touch all cuts present in the HPLs and 
2dHPLs on at least one of their corners, which is not the case for the (b,c,d)-type regions. For this reason, 
it is not possible to express the HPL and 2dHPL from (la) in terms of HPL and 2dHPL with same set of 
elements in the index vector in the (b,c,d)-type region. The case just considered of the continuation from 
(la)_ to (ld)_ is in this respect very fortunate, since the 'missing' cut in z = 1 (which is not touched by 
any boundary of (ld)_) translates into a new cut in s = —1, which requires only an extension of the set 
of indices of the HPL from (0,1) to (—1, 0, 1), while leaving the index set of the 2dHPL unaltered. The set 
(—1,0, 1) coincides with the original definition of the HPL and the numerical implementation jl^ is 
covering this set. 

It is clear that identities for the exchange and redefinition of arguments of the 2dHPL, which mix HPL 
and 2dHPL in a given region, owing to the presence of the value (—1) in the vector of the indices of the 
HPLs, are no longer applicable in (Id) _. An important consequence of this is furthermore that only a 



choice of variables of the form (5^) allows all functions in (ld)_ to be expressed in terms of HPLs and 
2dHPLs without an extension of the rational factors and the indices of the considered 2dHPLs. Hence, 
continuation to (lb,c)_ cannot be performed in two ways (as in the previous section), but only by using 
the second method, i.e. redefining the independent variables in the Euclidean region. 
In region (lb)_, (0 < z < 1, 1 — z < y < 1), we have 

< 0, Si2 < 0, < S23 < 0, -S23 > Sl3 > , (5.14) 

so that imaginary parts have to be assigned to S13 only. 

Continuation from (la)_ to (lb)_ is performed by first re-expressing the HPLs H(m;z) and 2dHPLs 
G{m{z);y) in (la)_ by HPLs H(m; x) and 2dHPLs G{rn{x);y), x = 1 — y — z, which is made according to 
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(3.6). The algorithm described in this section is then used to obtain the continuation to (lb)_ in terms of 
HPLs H(m;sib) and 2dHPLs G(TO(sib); rib), with 



r-ib 



^12 + ^23 

g2 



Slb 



£12 

g2 



(5.15) 



which fulfil < sib < 1, < rib < 1 — sib in (lb)-. For the reasons stated above, it is not possible to find 
a set of variables in (lb)_, which makes the (y <-> z) symmetry in this region explicit and still retains the 
same set of indices for HPLs and 2dHPLs. 

Finally, in region (lc)_, (0 < z < 1, — z < ?; < 0), we have 



< 0, Si3 < 0, < S23 < 0, -Si3 > Si2 > 



(5.16) 



so that imaginary parts have to be assigned to Si2 only. 

Continuation from (la)_ to (lc)_ follows similar lines by re-expressi ng t he HPLs H(to; z) and 2dHPLs 
G{rn{z);y) in (la)_ by HPLs H(m;y) and 2dHPLs G{rh{y); z), using (3.5). The algorithm described in 
this section is then used to obtain the continuation to (Ic)- in terms of HPLs H(m; sic) and 2dHPLs 
G(m(sic); ric), with 

Sl2 + Si3 Si3 

ric = = 1-2;, 



q^ 



Sic 



-y 



(5.17) 



which fulfil < Sic < 1, < ric < 1 — sic in (Ic)-. 

5.2 Continuation from (la) to (4d) and to (2d) ,(3d) 

In region (4d)_, (1 < j/ < oo, ~oo < z < — y), or {—oo < z < —1, 1 < y < — z), the invariants fulfil 

q^ < 0, si2 < < 0, si3 < < 0, S23 > -q^ > 0, 



(5.18) 



such that only S23 is to be assigned an imaginary part +ie. Equation (3.4) reads as above in Section ^.1: 

si3 S23 + 



^(312,313, S23 + ie) = (-g^)"^' 



S12 + si3 + S23 + «e S12 + si3 + S23 + le 



= {-qT'^iy + ie,z-ie) 



as, because of Eqs. ( 5.18 ) and by definition of y, z in Eqs. (2.5), 



Sl3 



Sl3 



S23 + «e S23 . 

„ + le ^ y + le , = — le ^ z - le . 

S12 + Si3 + S23 + le q-^ S12 + Si3 + S23 + «e q- 

To continue from (la)_ to (4d)_, we introduce new dimensionless variables, r^^ and s^d- 



rid 



si3 + S23 y + z 



Sid 



q 



1 



S23 Z S23 z 

and we express y and z in region (4d) in terms of rid a-nd 54^: 

1 

S4d 



1 - 

y = , z 



Sid 



(5.19) 
(5.20) 

(5.21) 
(5.22) 



As in all cases discussed before, rid-,Sid vary in the ranges < Sid < 1 and < r4d < 1 — Sid, i.e. the 
above parametrization maps region (4d)_ into region (la)_. In these variables, the generic expression for 
a master integral, Eq. (3.4), becomes 



«'(si2, si3, S23 + ie) = (-g^)"$ ' ^ ^^'^ 



1 



le, le 

Sid Sid 



(5.23) 
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HPLs and 2dHPLs at weight w = 1 thus become in (4d): 



H(0;z-ie) 
n{l;z-ie) 
GiO;y + ie) 
G{1- y + ie) 

G{1 - z + ie;y + ie) 
G(-z + ie-y + ie) 



log ( ze ) = - log (s4d) - 

V ■S4d / 

-H(0; S4d) - ITT 
-log (l + — +ie 
-H(-l;.S4d)+H(0;s4d) , 
log + 

V ■S4d 

G(l;r4d)-H(0;s4d) 
1 - Rd 



log 1 



S4d 



log 



1 - ?'4d ~ S4d 
S4d 



G(l - S4d; r4d) - H(l; S4d) - H(0; S4d) - in 



log 



l-''4d 



le + h «e 

^4d 



1 + ^ 



G(-S4d; f4d) + H(0; S4d) - H(-l; S4d) 
/ i-Tid , 1 

log 



ze te 

S4d 



= G(0;r4d). (5.24) 

The absence of imaginary parts from the last formula is non-trivial, and can be shown by explicitly inserting 
the definitions (5.21), taking account of the boundaries (5. IS) on the invariants. 

Continuation of the HPLs is made by combining the transformation formulae of Sections 5.1 (negation of 
argument) and 5.2 (inversion of argument) of |13|. Continuation of the 2dHPLs requires first the separation 
of all leftmost (— 2;)'s in the index vector. The remaining 2dHPLs are then continued according to 



G{a{z — ie), b{z — ie); y + ie) = 

= G{a{z — ie),b{z — ie); —z + ie) + 



y+ie 

y 



dy'g{a{z-ie),y')G{b{z-ie);y') 



G(a(z — ie), 6(z — ie); — z + ie) + / dy'g(a(z — ie), + ie)G(6(z — ie); y' + ie) 



(5.25) 



The term G(a(z — ie),b{z — ie);—z + ie) can be obtained by first evaluating G{a{z),b{z);—z), for < 
z < 1, according to the algorithm described in the Appendix of |^ and then continuing that result to 
z = — l/s4(j — ie, yielding an expression containing HPL H(m; S4d). For the second term in Eq. ( 5.25| ) one 
can introduce the integration variable r' = 1 — S4d2/', or y' = (1 — r')/s4d, so that expressing y, z in terms 
of r4d, S4d, Eqs. ( ^.22 ), the above equation becomes 



1 



G(a(z — ie), b{z — ie); y + ie) = G ( a ( ie ) , 6 ( ) ; H le 



'■4d 2 

dr' g I a 

Sid 



Sid 



Sid 

l-r' 

Sid 



1 



G b 



Sid 



1 

S4d 

l-r' 

S4d 



(5.26) 



where the expression for G(6(z); (1 — r\^)/ Sid) in terms of G(c(s4d); r'^^) and H(m; S4d) is again of lower 
weight and thus known. The expressions for the rational fractions are 



1 

! 

S4d 

1 

— \ 

Sid 



r{0;y') = -g(i;^4d), 
?(i;i/') = -g(i - S4d;<d) . 
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— g(l - z; y') = -g(-S4d; »-4d) : 

S4d 

—g{~z;y') - -g{0;r',^), (5.27) 

such that the r4j-integral in ( 3.25| ) yields a 2dHPL of argument rid- 
In region (2d)_, (1 < y < oo, 1 < z < oo), we have 

< 0, si3 < < 0, S23 <q^ <0, si2 > -q^ > 0, (5.28) 

so that imaginary parts have to be assigned to si2 only. For the continuation from (la)_ to (2d)_, we use 
again in (la)_ the x — 1 — y — zor si2 ^ Si3 interchange of Eq. (^^), so that its continuation to (2d)_ is 
given by 

*(S12 + S23) = (-9^)"$" ( ^ — , ^'^^t'' ^ ) = + ^ - • 

VS12 + Sl3 + S23 + «e Si2 + Sl3 + S23 + J 

(5.29) 

The algorithm already described in this section is then used to express $"(j/ + ie,x — ie) in the region (2d)_ 
in terms of HPLs H(m; S2d) and 2dHPLs G(TO(s2d); f'2d), with 

S12 + S13 x + y q^ 1 
?'2d = , S2d = = , (5.30) 

Sl2 X Sl2 X 

which fulfil < S2d < 1, < r2d < 1 - S2d in (2d). 

Finally, in region (3d)_, (— cx) < y < —1, 1 < z < — y) or (1 < z < oo, —oo < y < — z), the invariants 
are bound by 

q^ < 0, si2 < 9^ < 0, S23 < 9^ < 0, sig > -q^ > 0, (5.31) 

so that imaginary parts have to be assigned to S13 only. For the continuation from (la)_ to (3d)^, we use 
again in (la)_ the y ^ z or S13 <-> S23 interchange of Eq. ( |3.5| ), so that its continuation to (3d)_ is given 

by 

*(si2, S13, S23 + ^e) = (-g^)"*' f — -^f — , , ) = {-q^r^'{z + te,y- ze) . 

V S12 + Si3 + S23 + S12 + Si3 + S23 + «£/ 

(5.32) 

The algorithm described in this section is then used to express $'(z + ie,y — ie) in the region (3d)_ in 
terms of HPLs H(m; S3d) and 2dHPLs G(m(s3d); r^d), with 

si3 + S23 y + z q^ 1 

rsd = = , S3d = = — , (5.33) 

si3 y Si3 y 

which fulfil < S3d < 1, < r3d < 1 - S3d in (3d). 

5.3 Continuation from (la)_ to (4b) and to (2b)_, (3b)_ 

In region (4b) _, (l<z<oo,0<y + z< 1), the invariants fulfil 

q^ < 0, S12 < 0, S23 < g^ < Si3 < -S23 , (5.34) 



such that only S13 is to be assigned an imaginary part +ie. Equation (3.4) then reads 



^-(512, S13 + ^e, S23) = {-qT<i> [ '''Y' ^ ■ ' ^ — -- ) = i-q'r^iy - ^ + ■ 

\Si2 + si3 + S23 + S12 + si3 + S23 + lej 

(5.35) 

For the continuation from (la)_ to (4b) _ we introduce the dimensionless variables r^y, and S4h- 

S13 + S23 y + z S12 1-y-z 

r4b = = , S4b = — = , (5.36) 

S23 z S23 z 
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such that 

V = ; z = . (5.37) 

f4b + S4b ?'4b + S4b 

As above < S4b < 1 and < r4b < 1 — S4b, thus mapping region (4d) into region (la). In these variables, 
the generic expression for a master integral, Eq. (3.4), becomes 

^-(512, S13 + ie, S23) = ( ^ - ^ + iA ■ (5.38) 

V r4b + S4b r'4b + S4b / 



The continuation of HPLs and 2dHPLs at weight w = 1 reads 

( 



H(0;z + ie) = log ( \- ie] = - log ( ^ ^^'^ ) + log S4b 

r4b + S4b / V S4b 



= -G(-S4b;r-4b) - H(0;s4b) 
11(1; 2: + le) = — log 1 le \ = — log 



?'4b + S4b / V ''4b + S4b 

= -G(l - S4b; 7'4b) + G(-S4b; r4b) + H(0; S4b) + H(l; S4b) + in 
G(0; y-ie) = log ( — - — ^ ie) = log(l - r4b) - log(r4b + S4b) - in 

\ ''4b + S4b / 

= G(l; r4b) - G(-S4b; ?'4b) - H(0; S4b) - in 
Gil■y-^e) = log fl + ^-111^ + log ^ ' + -^^-^ 



f 4b + S4b / V ^ib + S4b 

= -G(-S4b; r4b) - H(0; S4b) + H(-l; S4b) 

/l + ^-^+ie i i€\ I SA^ 

G{1- z-ie-y-ie) = log '^^i^ j ^-^^^ = log ' ^'^ 



1 - / V 1 - 7'4b - S4b 

■r4b + S4b 



= -G(l - S4b; r4b) + H(0; S4b) + H(l; S4b) + in 

(~-^=^-ie+ } +ie\ 

G{-z -ie-y~ le) = log ''^^+^^^ , ""-^^ = log(r4b) 

\ T h / 

- G(0;r4b). (5.39) 

At variance with all cases discussed before, one observes here that the 2dHPLs G(6(s4b); r4b) appear not 
only in the continuation of the 2dHPLs G(6(z); y), but also in the continuation of the HPLs H(&; z). This 
feature is due to the fact that r4b appears in the expressions for both y and z ( ^.37 ), while it appeared only 
in the expression for y in all cases discussed previously. As a consequence, the continuations of G{b{z);y) 
and H(6; z) of weights w > 1 are more intertwined than in the cases discussed in previous sections, and no 
simple formulae for them can be given. Instead, these continuations have to be carried out in an algorithmic 
procedure, which is explained below. 

Continuation of the HPLs is made using 



H(o,6;z + ie) = / dz' f(a, z')H(6; z') 
Jo 



dz' f(a, z')H(6: z') + f ^ dz' f(a, z')H(6; z') (5.40) 

= 4b 

H ( a, fe; — + V . , '^'^^^ .. f ( a: — ^ + le ] U \ b 



S4h ' "V Jo ('■4b + 'S4b)^ A"' ''4b + *4b ' " V' ''4b + *4b 

where the new integration variable r^^ was introduced with the substitution z' — l/(?'4b + Sih) + One 
finds for the rational fractions: 

^ M 0; , „ ) = g(-S4b,r4b) 



(''4b + S4b)^ V ''4b '^4b 
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K,h^'{'-^7r^) = g(i-4.,.;.)-g(-.4.,.l.). (5.41) 



The HPLs H(6; l/{r'^i, + S4b) + *e) of Eq. (5.40) are of lower weight, and therefore known when proceeding 
bottom up starting from weight w — 1. They can be expressed as a hnear combination of HPLs H(c; s^h) 
and 2dHPLs G{d{s4h),r'^^)- As a consequence, the above integral yields 2dHPLs G(e(s4b), ?'4b). Finally, 
the boundary term H(a, b; l/s4b + ?e) is evaluated using the inversion formula of Section 5.3 of p^ , yielding 
HPLs H(c;s4b). 

For obtaining the continuation of the 2dHPLs, write 

G{c{z + ie);y-ie)^G(c( ^ + ie);- ^ ~ '^^'^ - if) , (5.42) 













; — — -ie) 


\ \Sih 




S4b / 



7'4b + S4b / ''4b + S4b 

and observe that the r.h.s., considered as a function of r4b, is equal to its value at r4b = plus the integral 
of its derivative from to r4b, i.e. 

G(c(z + ie);y- 

\ \ «4b 

dr4bT^Gfcf— ^ +ie]--4^^-ie] . (5.43) 

The r^j^-derivative in the above formula acts both on the argument of the 2dHPL and on the index 
vector. Writing out the 2dHPL in its multiple integral representation, this derivative can be carried out, 
yielding at most the squares of inverse rational factors. Using partial fractioning and integration by parts, 
the result of this differentiation can be rewritten as a linear combination of integral representations of 
2dHPLs G((i(s4b), r4b); the algebraic simplifications occurring in working out the arguments of the f actors 
;(a,r4jj), appearing in the r^i^-derivatives are similar to those encountered in the r.h.s. of Eq. ( 5.39| ). The 



boundary term G(c(l/s4b + ie); — S4b ^ *e) is obtained by first working out G(c(z); — z) for < z < 1 by 
the standard procedures described in the appendix of then continuing the result to z = — l/s4b; it is 
expressed in terms of HPLs H((i; S4b). 

In region (2b)_, (0 < z < 1, — oo < y < —z), 

< 0, S23 < 0, si2 < g^ < si3 < -si2 , (5.44) 
so that imaginary parts have to be assigned to S13 only. Continuation from (la)_ to (2b)_ uses again the 



a: = l— y — zor S12 ^ S13 interchange of Eq. ( |3.6| ), as in the continuation from (la)_ to (2d)_, Eq. (5.29). 
Subsequently, the algorithm described in this section is then used to obtain the continuation to (2b) _ in 
terms of HPLs H(m; S2b) and 2dHPLs G(m(s2b); ?'2b), with 

S12 + S13 x + y S23 l-x~y 
r2b = = , S2b = — = , (5.45) 

S12 X S12 X 

which fulfil < S2b < 1, < r2b < 1 - S2b in (2b). 

Finally, in (3b)_, {1 < y < 00, —y < z < —y + 1), the invariants are bound by 

< 0, S12 < 0, si3 < g^ < S23 < -S13 , (5.46) 

so that imaginary parts have to be assigned to S23 only. Continuation from (la)_ to (3b)_ employs again 



the y ^ z interchange of (3.5), as in the continuation from (la)_ to (3d)_, Eq. (5.32). The algorithm 
described in this section is then used to obtain the continuation to (3b)_ in terms of HPLs H(to; S3b) and 
2dHPLs G(TO(s3b);r3b), with 

S13 + S23 y + z S12 l-y-z 

r-ih = = , S3b = — = , (5.47) 

si3 y si3 y 

which fulfil < ssb < 1, < rsb < 1 — S3b in (3b). 
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5.4 Continuation from (la)_ to (3c)_ and to (4c)_ and (2c)_ 

The last regions required for the kinematics of DIS are (3c)_ and (4c)_, (2c)_, related to it by crossings. 
In region (3c)_, (!<?/< oo, < z < 1) the invariants are bound by 

< 0, si2 > 0, si3 < < 0, < S23 < , (5.48) 



such that only Si2 is to be assigned an imaginary part +ie. Eq. (3.4) then reads 

*(si2 + le, si3, 523) = ( V)"*!* ( — , — —) = i-q"THy + ^ + *e). 

\Sl2 + Si3 + S23 + «e Si2 + Si3 + S23 + 

(5.49) 

For the continuation from (la)_ to (3c) _ we introduce the dimensionless variables are r^c and s^c' 

Sl2 + S23 y - 1 S23 z 

rsc = = , S3c = — = - , (5.50) 

si3 y si3 y 

such that 

z = -^. (5.51) 

1 - 7'3c 1 - f 3c 

As above, < ssc < 1 and < r3c < 1 — ssc, thus mapping region (3c) into region (la). In these variables, 
the generic expression for a master integral, Eq. (3.4), becomes 

■^{si2 + le, si3, S23) - (-?')"$ (—^ + ie, + ze) . (5.52) 



, 1 - 7-3C 1 - r3, 

The HPLs and 2dIIPLs at weight w — 1 are continued as: 
H(0;z + ^e) - fog f-^^ 

V 1 - 7-30 

= H(0;s3c)-G(l;r3c) 
H(l;z + ie) = -logfl 



G{Q-y + ie) = log 



1 - 7'3c 

-G(l - S3c; rgj + G(l; rgc) + H(l; S3c) 
1 



1 - n3c 

= -G(l;r3c) 
G(l;y + «e) = log ( 1 - — ) = log 



1 - ?'3c / V 1 - ''3c 

- G(0;r3c)-G(l;r3c)-i^ 
G{l-z-^e■y + ^e) = log f ^ ~ ^ 'J' ~ ' ^ fog f /^^ + ) - ^^ 

= G(-S3c; r3c) - G(l - S3c; r3c) + 11(0; s^^) + H(l; s^^) - ^ , 



^3c 



G(-z -ie-v + le) = log I ^''^^^ ^''^^ I = log 



1 + S3c 



\ l-r3c + / V 53c 

= H(-l;s3c)-H(0;s3c) . (5.53) 

As in the continuation to (4b)_, one observes here that 2dIIPLs G(6(s3c); r^c) appear both in the continu- 
ation of the 2dHPLs G{b{z);y) and of the HPLs 11(5; z). The continuation of the higher-weight functions 
also follows similar lines as for (4b) _. 
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As in Eq. (|5.40|) , continuation of the HPLs is made using 



H(a, b;z + ie) ^ H (a, b; sgc + ie) + [ dz'f(a, z')Ti{b; z') (5.54) 



H( 



One finds for the rational fractions: 



' 3c 

= g(l;r^J-g(l-S3,;r;,J . (5.55) 

\ 3c / \ 3c / 

The HPLs H(6; s^c/ (1 — ^3c) + ^-re of lower weight, and therefore known. They can be expressed as linear 

combination of HPLs H(c;s3c) and 2dHPLs G((i(s3c), rg^.). As a consequence, the above integral yields 
2dHPLs G(e(s3c),r3c). 

To continue the 2dHPLs, following Eqs. (|!4|) and (|^) write 

G{b{z + ie)-y + ie) = G (b ( + ie] ; — + leA (5.56) 



, 1 - f'3c / 1 - 7'3c 

and observe that the r.h.s., considered as a function of rgc, is equal to its value at rgc — plus the integral 
of its derivative from to rgc, i.e. 



G{b{z + ie):y + ie) = G (ssc + ie) ; 1 + ie) 



d^ScX^Glfel-^ + zeV-i— +ze) . (5.57) 



d^^Sc V V 1 - 7-3C / ' 1 - ?-3c 

As in the previous section, the rg^-derivative in the above formula acts both on the argument of the 
2dHPL and on the index vector. Writing out the 2dHPL in its multiple-integral representation, this deriva- 
tive can be carried out, yielding at most the squares of inverse rational factors. Using partial fractioning 
and integration by parts, the result of this differentiation can be rewritten as a linear combination of inte- 
gral representations of 2dHPLs G(c(s3c), rgc). The boundary term G(6(s3c -I- ie); 1 -I- ie) is obtained using 
standard procedures, as described in the appendix of yielding H(c;sgc). 

In region (4c)_, (0<?;<1, 1 < z < oo), 

< 0, si2 > 0, S23 < < 0, g2 < sig < , (5.58) 

so that imaginary parts have to be assigned to si2 only. Continuation from (la)_ to (4c)_ uses again the 
z ^ y interchange of ( [4.3l| ). Subsequently, the algorithm described in this section is then used to obtain 
the continuation to (4c)_ in terms of HPLs H(m; S4c) and 2dHPLs G(to(s4c); r4c), with 

512 + Si3 2-1 Si3 y 

= = , S4c = — = - , (5.59) 

S23 Z S23 Z 

which fulfil < S4c < 1, < r4c < 1 — S4c in (4c). 

Finally, in region (2c)_, {0 < y < 1, —oo < z < —y), the conditions on the invariants are 

q^ < 0, S23 > 0, si2 <q^ <0, q^ <si3<0 , (5.60) 

so that imaginary parts have to be assigned to S23 only. Continuation from (la)_ to (4c) _ employs as 
well the y ^ z interchange of (4.31), followed by the z = 1 — x — y replacement of ( 1.2l| ), which maps 
(4c) into (l<a:<c», Q < y < \)\ note that this procedure, which involves both the interchange and 
the replacement of the arguments, differs from all crossings discussed up to here. The algorithm described 
in this section is then used to obtain the continuation to (2c) _ in terms of HPLs H(m; S2c) and 2dHPLs 
G(m(s2c);'^2c), with 

513 + S23 a;-l si3 y 

r2c = = , S2c = — = - , (5.61) 

S12 X S12 X 

which fulfil < S2c < 1, < r2c < 1 — S2c in (2c). 
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Region 


Variab 
y-type 


es 

z-type 


Procedure 




7,1 - ^ 

yia - ^ 


2la — T 




[ID) 


rib 


q 


l,i/la, 2:iaj > (.i/la, 1 i/la ^laj ^ I'lbiSlbj 


[IC) 


' Ic — 

q 


Ale — 

q 


l,i/la,2la; (^2la,J/la; * {'Ic^lc) 




„ ^ _ Sl2 + S23 
' Id 2 

q 


Sid - ^2 


( 11-1 1 1 V 1 T^T 1 Qt 1 1 

l,i/la,^la; * Vld5Sld; 




U2a - si2 


q^ 

"2a - 


\ 1 ti '/ ^ \ - ^ [ill \ - 1 i■^ - - 1 - J. (ill-. "7 1^-. 1 

l,i/lai2ia^ > \yia, i- yia. ^la) ' {"■2a,V2a) 


(ZD) 


_ Sl2 + Sl3 


S2b - 


l,i/lai2ia^ > \yia, ^ yia ^la) ' l,'2b;S2b; 


[ZC) 


^_ _ Sl3 + S23 


St - £13 
^2c - si2 


l,i/lai2laj 1^ (^2ia,(/la; ' {^la, ^ yia ^la) ' \i2c;02c) 


(za) 


r J - ''^12 + Sl3 


2 

S2d - si2 


l,i/la, ^la; * li/la, 1 2/la ^laj ^ I' 2d, S2d ; 




M3a - 


g2 

W3a - 


lilt "yt 1 i. ( "y t lit 1 i. (ill-. Tin 1 

l,i/la,2ia^ > \Zia,yia) * \^J-Zai^JZa) 


(6D) 


r,u — •^13 + *23 
'^3b - si3 


Sl2 

S3b - 


lilt 'y t 1 ( ^ t lit 1 ( 1^ f-\1 1 

l,J/la,2la; > (,21a, J/la; > (,'3b,S3b; 




_ Si2 + S23 
'3c - si3 


- £23 
«3c - 


lilt "y t \ i ( ,-, o „ 1 

l,i/la,2laj * (,'3c,S3cj 


[M) 


_ Si3 + S23 


2 

«3d - 


lilt "y t \ i ( "y t lit 1 i ( T'r^ 1 e n 1 1 

(yiajZla) > (,21a, ylaj * (,'3d,S3dj 




UA - 

M4a - S23 


q^ 

W4a - 


(,i/la,21a; ' V'^4a, t74a; 


(4b) 


_ Si3 + S23 
'^4b S23 


- S23 


(yia,2ia) (r4b,S4b) 


(4c) 


^ _ Si2 + Sl3 
'^"c - S23 


- S23 


(yia,2la) {Zla,yia) (»"4c,S4c) 


(4d) 


J,, , _ Sl3 + S23 
'^^d - S23 




(yia, Zia) {rid, Sid) 



Table 1: Variables used for 2dHPLs and HPLs in each kinematic region, and variable transformations 
applied in the analytic continuation from (la) to each region. 

6 Conclusions 

In this paper, we have described the analytic continuation of two-loop four-point functions with massless 
internal propagators and one off-shell external leg from the Euclidean region to all Minkowskian regions 
of physical interest. While the continuation to the region relevant to 1 — > 3 kinematics (as in e+e" 
3 jets) amounts to simply continuing an overall factor ^, the continuations for the 2 — > 2 scattering 
processes {V -\- Ij and DIS-(2 -I- l)j production) is more involved. In particular, since all two-loop four- 
point functions are expressed in terms of 2dHPLs whose arguments lie in general outside the analyticity 
range, for which numerical routines are available jl^ , we had to find appropriate variable substitutions to 
map each kinematical region into the analyticity range of the 2dHPLs. These variable transformations are 
summarized in Table |^. 

Each continuation of the two-loop four-point functions is performed by first identifying the variables 
crossing a kinematical cut. Subsequently, the basis functions of the two-loop four-point functions (HPLs 
and 2dIIPLs) are continued by first continuing the weight w — 1 functions (which are just logarithms 
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with a known and well-defined continuation), the functions of higher weight are then constructed using 
the product algebra and the integral representations of the HPLs and 2dHPLs. It turns out, using this 
approach, that all imaginary parts arising in the analytic continuation are made explicit, and their signs 
are fixed from the continuation of the w = 1 functions, which are all listed in the appropriate sections. 

Using the continuation formulae derived in this paper, one can use the results obtained for the two- 
loop QCD corrections to the e+e^ 3 jets matrix element and hclicity amplitudes to derive the 
corresponding corrections to the matrix elements for vector-boson-plus-jet production at hadron colliders 
and to deep inelastic two- plus-one-jet production [po||. This work is currently in progress. 



A Harmonic polylogarithms 

The generalized polylogarithms Sn,p{x) of Nielsen [p^ turn out to be insufficient for the computation 
of multi-scale integrals beyond one loop. To overcome this limitation, one has to extend generalized 
polylogarithms to harmonic polylogarithms 

Harmonic polylogarithms are obtained by the repeated integration of rational factors. If these rational 
factors contain, besides the integration variable, only constants, the resulting functions are one-dimensional 
harmonic polylogarithms (or simply harmonic polylogarithms, HPLs) If the rational factors depend 

on a further variable, one obtains two-dimensional harmonic polylogarithms (2dHPLs) In the 

following, we recall the definition of both classes of functions, and summarize their properties. 

A.l One-dimensional harmonic polylogarithms 

The HPLs, introduced in [p^ , are one- variable functions H(a; a;) depending, besides the argument x, on a 
set of indices, grouped for convenience into the vector a, whose components can take one of the three values 
(1,0,-1) and whose number is the weight w of the HPL. More explicitly, the three HPLs with w — 1 are 
defined as 



H(l;x) = / ^ = -ln(l-x) 

Jo ^ — X 

H(0; x) = Inx , 



dr' 

H(-l;x) = / ^--^=ln(l+x) ; (A.l) 



their derivatives can be written as 

d 



H(a; x) = f(a; x) , a = 1,0,-1, (A.2) 
da; 

where the 3 rational fractions f(a; x) are given by 

f(0;x) - - , 

X 

f(-l;x) = (A.3) 
I + X 

For weight w larger than 1, write a — (a, 6), where a is the leftmost component of a and h stands for the 

vector of the remaining {w — 1) components. The harmonic polylogarithms of weight w are then defined 
as follows; if all the w components of a take the value 0, a is said to take the value 0^ and 

H(0^;a:) = ^In'^a; , (A.4) 

while, if a 7^ 0„, 



H(a; x) 



[ dx' i{a;x') R{b;x') . (A.5) 
Jo 
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In any case the derivatives can be written in the compact form 



4^H(a; x) = f(a; a;)H(6; x) , (A.6) 
dec 

where, again, a is the leftmost component of a and b st and s for the remaining [w — 1) components. 

It is immediate to see, from the very definition Eq. (A^), that there are 3"" HPLs of weight w, and that 
they are hnearly independent. The HPLs are generalizations of Nielsen's polylogarithms |Q. The function 
Sn,p{x), in Nielsen's notation, is equal to the HPL whose first n indices are all equal to and the remaining 
p indices all equal to 1: 

) =H(0„,lp;x) ; (A.7) 
in particular the Euler polylogarithms Li„(x) = S„-i,i(a;) correspond to 

Li„(a;) =H(0„_i,l;x) . (A.8) 

As shown in , the product of two HPLs of a same argument x and weights p, q can be expressed as 
a combination of HPLs of that argument and weight r = p + q, according to the product identity 

H(p;x)H(9;x) = ^ H(f; x) , (A.9) 

r— pUlg 

where p , (f stand for the p and q components of the indices of the two HPLs, while pb) q represents all 
mergers of p and q into the vector r with r components, in which the relative orders of the elements of p 
and q are preserved. 

The explicit formulae relevant up to weight 4 are 

H(a; a;) H(6; x) = H(a, b; x) + H(6, a; x) , (A.IO) 



H(a; x) H(6, c; x) 
H(a; x) H(6, c, d; x) 



H(a, 6, c; x) + H(6, a, c; x) + H(6, c, a; x) , 

H(a, 6, c, d; x) + H(6, a, c, d; x) + H(6, c, a, d; x) + H(6, c, d, a; x) , (A. 11) 



and 



H(a, 6; x) H(c, d; x) = H(a, 6, c, d; x) + H(a, c, 6, d; x) + H(a, c, d, 6; x) 
+ H(c, a, 6, d; x) + H(c, a, d, b; x) + H(c, d, a, 6; x) 



(A.12) 



where a,b,c,d are indices taking any of the values (1,0,-1). The formulae can be easily verified, one at 
a time, by observing that they are true at some specific point (such as x = 0, where all the HPLs vanish 
except in the otherwise trivial c ase i n which all the indices are equal to 0) , then taking the x-derivatives of 
the two sides according to Eq. ( A.6 ) and checking that they are equal (using when needed the previously 

established lower- weight formulae). 

Another class of identities is obtained by integrating (A.4) by parts. These integration-by-parts (IBP) 
identities read: 



H(mi 



,TOq;x) 



H(mi;x)H(rn2, ■ • ■ ,mq;x) - H(m2, mi; x)H(m3, . . . ,mq\x) 



'97 ■ 



, mi;x) 



(A.13) 



These identities are not fully linearly independent of the product identities. 

A numerical implementation of the HPLs up to weight ui = 4 is available [O] . 
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A. 2 Two-dimensional harmonic polylogarithms 

The 2dHPLs family is obtained by the repeated integration, in the variable y, of rational factors chosen, in 
any order, from the set 1/y, \/{y ~ 1), 1/(2/ + z — 1)^ l/(y + -2), where z is another independent variable 
(hence the 'two-dimensional' in the name). In full generality, let us define the rational factor as 

g{a-y)^-^, (A.14) 

y- a 

where a is the index, which can depend on z, a = a(z); the rational factors that we consider for the 2dHPLs 
then are 



g(0;2/) 
g(i;2/) 

;(1 - z]y) 



y + z-1 



gi-z;y) = . (A.15) 

y + z 

With the above definitions, the index takes one of the values 0, 1, (1 — z) and (—z). 

Correspondingly, the 2dHPLs at weight w = 1 (i.e. depending, besides the variable y, on a single further 
argument, or index) are defined to be 

G(0;y) = Iny, 
G(l;y) = In(l-y), 

y 



G{1 -z;y) = In 1- 



1 



G(-z;y) = ln(l + ^V (A.16) 



z > 

The 2dHPLs of weight w larger than 1 depend on a set of w indices, which can be grouped into a w- 
dimensional vector of indices a. By writing the vector as a = (a, 6), where o is the leftmost component of a 
and h stands for the vector of the remaining (w — 1) components, the 2dHPLs are then defined as follows: 
if all the w components of a take the value 0, a is written as 0^ and 

G{^^-y) = -yy , (A.17) 

while, if a 7^ Ouj, 

G(a;y)= [' dy' g{a;y') G{b;y') . (A.18) 
Jo 

In any case the derivatives can be written in the compact form 

^G{a;y)^gia;y)G{b;y) , (A.19) 
dy 

where, again, a is the leftmost component of a and b stands for the remaining [w — 1) components. 

It should be observed that the notation for the 2dHPLs employed here is the notation of [Q, which 
is different from the original definition proposed in H]. Detailed conversion rules between the different 
notations, as well as relations to similar functions in the mathematical literature (hyperlogarithms and 
multiple polylogarithms) can be found in the appendix of . 

Algebra and reduction equations of the 2dHPLs are the same as for the ordinary HPLs. The product of 
two 2dHPLs of a same argument y and weights p, q can be expressed as a combination of 2dHPLs of that 
argument and weight r = p + q, according to the product identity 

G{p;x)Giq;x) = ^G(r;x), (A.20) 

r—p[+iq 
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where p, q stand for the p and q components of the indices of the two 2dHPLs, while p^ q represents all 
possible mergers of p and (f into the vector r with r components, in which the relative orders of the elements 
of p an d q a re preserved. The explicit product identities up to weight w = A are identical to those for the 
HPLs ( |A.10D -( |a!T2|) , with all H replaced by G. 
The integration-by-parts identities read: 



G(mi, . . . ,mq; x) 

A numerical implementation of the 2dHPLs up to weight w = 4 is available |1 



G{mi;x)G{m2, . . . , ruq] x) - G(m2, mi; a;)G(m3, . . . , m,; x) 
+ ... + {~iy+^G{mq, ...,mi;x) . (A.21) 
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